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Abstract.  This  paper  describes  the  methodology  used  by  the  General  Rotorcraft 
Aeromechanical  Stability  Program  (GRASP)  to  model  the  characteristics  of  the  flow 
through  a  helicopter  rotor  in  hovering  or  axial  flight.  Since  the  induced  flow  plays  a  sig¬ 
nificant  role  in  determining  the  aeroelastic  properties  of  rotorcraft,  the  computation  of 
the  induced  flow  is  an  important  aspect  of  the  program.  Because  of  the  combined  finite- 
element /multibody  methodology  used  as  the  basis  for  GRASP,  the  implementation  of 
induced  velocity  calculations  presented  an  unusual  challenge  to  the  developers.  To 
preserve  the  modeling  flexibility  and  generality  of  the  code,  it  was  necessary  to  depart 
from  the  traditional  methods  of  computing  the  induced  velocity.  This  is  accomplished 
by  calculating  the  actuator  disc  contributions  to  the  rotor  loads  in  a  separate  element 
called  the  air  mass  element,  and  then  performing  the  calculations  of  the  aerodynamic 
forces  on  individual  blade  elements  within  the  aeroelastic  beam  element.  ( ■  (2  ^  w  6 1  '  3 

Keywords.  Aeroelasticity;  Finite  Elements;  Helicopters;  Induced  Velocity;  notary 
Wings;. 
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INTRODUCTION 

In  September  1980,  work  began  on  developing  the 
General  Rotorcraft  Aeromechanical  Stability  Program 
(GRASP).  While  numerous  analyses  (Ormiston  and 
Hodges,  1972;  Friedmann,  1973;  Hodges,  1976,  1979; 
Warmbrodt  and  Friedmann,  1979;  Friedmann  and 
Straub,  1980;  Davis  et  of.,  1974;  Bielawa,  1976;  John¬ 
son,  1977,  1980;  Sivaneri  and  Chopra,  1982)  are  avail¬ 
able  to  perform  aeroelastic  analyses  for  rotorcraft,  all  of 
them  are  subject  to  major  limitations  (Johnson,  1986) 
in  generality,  flexibility,  or  theoretical  consistency.  The 
purpose  for  which  GRASP  has  been  developed  is  to 
provide  a  tool  with  enhanced  capabilities  that  can  be 
used  to  perform  aeroelastic  calculations  for  helicopters 
in  hover  and  axial  flight. 

The  implementation  of  the  hybrid  finite-element/multi- 
body  methodology  (Hodges  et  of.,  1987a,  1987bb)  in 
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GRASP  allows  a  structure  to  be  modeled  as  a  collec¬ 
tion  of  rigid  bodies  and  flexible  elements  that  can  be 
connected  in  a  completely  arbitrary  manner.  While  this 
methodology  presents  the  analyst  with  a  great  deal  of 
generality  and  flexibility  in  structural  modeling,  it  also 
presents  the  developer  with  some  challenges  in  imple¬ 
menting  an  appropriate  representation  of  the  helicopter 
flow  field.  Since  the  treatment  of  the  flow  around  and 
through  the  rotor  disk  is  an  important  part  of  any  aero¬ 
elastic  analysis  of  rotorcraft,  it  is  vital  that  the  induced 
velocity  be  calculated  in  a  consistent  manner. 

In  this  paper,  current  methods  used  to  calculate  the 
inflow  will  first  be  described.  Then,  the  methodology 
used  in  GRASP  will  be  discussed  and  the  differences 
with  the  more  traditional  methods  highlighted.  Finally, 
the  theoretical  basis  of  the  approach  implemented  in 
GRASP  will  be  outlined. 

METHODOLOGY 

Current  aeroelastic  stability  analyses  for  helicopters  use 
a  variety  of  methods  to  calculate  the  steady-state  and 


dynamic  induced  inflow.  These  range  from  simple,  lin¬ 
ear  models  for  uniform  inflow  to  sophisticated,  nonuni¬ 
form  inflow  models  using  free-wake  analyses.  While 
not  breaking  any  new  ground  with  respect  to  devel¬ 
oping  new  models,  GRASP  does  take  a  different  ap¬ 
proach  with  regard  to  its  calculations  of  induced  veloc¬ 
ity.  Therefore,  before  describing  the  methodology  used 
in  GRASP,  and  its  rationale,  it  will  be  instructive  to 
look  at  some  representative  examples  of  the  approaches 
taken  in  current  analyses. 

Traditional  Methods 

In  many  analyses  (Ormiston  and  Hodges,  1972;  Fried¬ 
mann,  1973;  Hodges,  1976,  1979;  Warmbrodt  and  Fried¬ 
mann,  1979;  Friedmann  and  Straub,  1980),  the  steady- 
state  induced  velocity  is  calculated  from  a  single,  linear, 
closed-form  expression  that  combines  both  momentum 
and  blade-element  contributions  to  the  rotor  forces  and 
moments.  This  expression  is  a  function  of  the  rotor 
collective  pitch  angle  (usually  at  the  three-quarter  ro¬ 
tor  radius).  Assuming  uniform  inflow,  one  takes  the 
induced  inflow  velocity  over  the  entire  rotor  disk  to  be 
constant  with  the  same  value  as  the  theoretical  value  at 
the  three-quarter  rotor  radius.  Alternatively,  one  could 
assume  that  the  inflow  angle,  which  is  the  inflow  veloc¬ 
ity  divided  by  the  local  blade  speed,  is  constant  over 
the  rotor  radius  with  the  same  value  as  the  theoretical 
value  at  the  three-quarter  rotor  radius. 

Another  approach  to  calculating  the  steady-state  inflow 
is  demonstrated  by  Davis  et  al.  (1974),  Bielawa  (1976), 
Johnson  (1977),  and  Sivaneri  and  Chopra  (1982).  In 
this  method  the  inflow  velocity  is  calculated  as  a  func¬ 
tion  of  the  thrust  coefficient,  which  is  usually  given. 
However,  the  blade  pitch  angle  required  to  produce  the 
desired  thrust  is  also  a  function  of  the  inflow  velocity. 
Thus,  the  computation  of  the  induced  velocity  is  nonlin¬ 
ear,  and  requires  an  iterative  solution.  The  distribution 
of  induced  velocity  over  the  rotor  disk  is  then  either  as¬ 
sumed  to  be  uniform,  or  specified  by  a  set  of  assumed 
functions  such  as  the  Glauert  induced  velocity  (Bielawa, 
1976). 

The  method  used  in  CAMRAD  (Johnson,  1980)  to  cal¬ 
culate  the  induced  velocity  is  more  sophisticated  than 
any  of  the  preceding  analyses.  CAMRAD  can  use  any 
of  three  methods  to  determine  the  induced  flow.  First, 
as  above,  a  uniform  inflow  distribution  is  computed  as 
nonlinear  function  of  the  thrust  coefficient.  Then,  if  de¬ 
sired,  a  nonlinear  distribution  can  be  determined  from 
a  prescribed-wake  analysis,  using  the  uniform  inflow  as 
an  initial  guess.  If  further  refinement  is  needed,  a  free- 
wake  analysis  is  performed  using  the  prescribed-wake 
solution  as  the  initial  guess. 

Of  the  analyses  just  discussed,  only  a  few  (Johnson, 
1977,  1980)  consider  the  effects  of  inflow  dynamics.  Ba¬ 
sic  to  this  type  of  dynamic  inflow  analysis  is  the  assump¬ 
tion  that  total  forces  on  the  rotor  vary  slowly  erough 
that  actuator  disk  theory  is  applicable  to  perturba¬ 
tion  velocities.  Comparisons  with  experimental  results 
(Johnson,  1986)  have  shown  that  dynamic  inflow  can 
have  a  significant  effect  on  aeroelastic  phenomena. 

A  common  feature  of  all  of  the  analyses  discussed  earlier 
is  that  the  calculation  of  the  steady-state  inflow  veloc¬ 


ities  is  performed  separately  from  the  main  calculation 
of  the  steady-state  deformation  of  the  structure.  Thus, 
the  inflow  generalized  coordinates  are  not  included  in 
the  state  vector  with  the  structural  degrees  of  freedom. 
However,  when  inflow  dynamics  are  included  in  the  dy¬ 
namic  problem  (Johnson,  1977,  1980),  those  generalized 
coordinates  are  (and  must  be)  included  in  the  state  vec¬ 
tor.  Although  this  is  a  somewhat  inconsistent  treatment 
of  the  steady-state  and  dynamic  inflow  generalized  co¬ 
ordinates,  it  does'  not  result  in  any  significant  analytical 
problems.  This  is  a  result  of  the  coupling  between  the 
steady-state  inflow  and  the  structural  deformations  be¬ 
ing  very  weak. 

To  separate  the  steady-state  induced  velocity  calcula¬ 
tions  from  the  structural  calculations,  contributions  to 
the  rotor  loads  from  individual  blade  elements  must  be 
calculated  at  the  same  time  as  the  flow  field  contribu¬ 
tions.  Since  the  exact  geometry  of  each  blade  may  not 
be  known,  it  is  necessary  to  assume  a  relationship  be¬ 
tween  one  or  more  blade  parameters  and  the  forces  on 
the  rotor.  For  example,  in  Sivaneri  and  Chopra  (1982) 
the  induced  velocity  is  calculated  as  a  function  of  thrust, 
and  thrust  is  a  function  of  the  blade  pitch  angle  at  the 
three-quarter  radius  and  the  induced  velocity.  It  is  ap¬ 
parent  in  this  case  that  there  is  some  implied  relation¬ 
ship  between  the  blade  pitch  angle  and  the  blade  geom¬ 
etry  and  section  aerodynamics. 

GRASP  Methodology 

The  axisymmetric  flow  field  for  a  helicopter  in  hover 
or  axial  flight  is  represented  in  GRASP  by  an  element 
called  the  air  mass  element.  The  inflow  generalized  co¬ 
ordinates  associated  with  this  element  are  then  included 
in  the  steady-state  model  state  vector  as  well  as  in  the 
dynamic  model  state  vector.  This  means  the  steady- 
state  inflow  velocity  is  calculated  in  parallel  with  the 
structural  deformations,  and  that  those  velocities  are 
fully  coupled  to  the  deformed  state  of  the  rotor  blades. 
The  inflow  generalized  coordinates  are  introduced  into 
the  model  in  a  manner  similar  to  that  used  to  introduce 
structural  degrees  of  freedom.  That  is,  an  air  node  is 
introduced  to  represent  the  flow  field  at  a  point  on  the 
axisymmetric  axis  of  that  flow  field. 

One  of  the  difficulties  that  arises  from  integrating  the 
flow-field  model  into  the  structural  model  is  related  to 
the  specification  of  the  motion  of  the  flow  field  rela¬ 
tive  to  the  structure.  It  is  known  that  if  the  rotor  disk 
undergoes  large  deformations,  the  flow  field  will  also  un¬ 
dergo  changes  that  follow  and  lag  the  disk  deformations. 
This  occurs  because  the  flow  field  is  not  physically  at¬ 
tached  to  the  rotor,  but  is  highly  dependent  (to  say  the 
least)  on  its  location.  However,  since  such  large  motions 
would  result  in  periodic  forces  and  moments  for  which 
GRASP  will  not  account,  that  situation  may  be  ignored. 
There  are,  therefore,  two  possible  implementations  that 
may  be  used.  The  first  associates  the  flow  field  with 
an  inertially  fixed  frame  of  reference.  In  this  case,  it  is 
understood  that  large  motions  of  the  rotor  (which  are 
not  allowed)  will  have  no  effect  on  the  location  of  the 
flow  field.  The  other  option  is  to  attach  the  flow  field  to 
the  structure  with  the  understanding  that  large  defor¬ 
mations  will  violate  the  flow-field  model  assumptions. 


For  GRASP,  the  former  option  was  chosen  since  it  is 
closer  to  the  actual  physics  of  the  phemonenon. 

Another  difficulty  with  integrating  the  air  mass  element 
into  the  structural  model  arises  because  of  the  mul¬ 
tilevel  substructuring  capabilities,  which  enhance  the 
flexibility  and  generality  of  GRASP  in  modeling  com¬ 
plex  structures.  One  of  the  concepts  fundamental  to  the 
use  of  multilevel  substructuring  is  that  no  substructure 
is  required  to  have  any  specific  knowledge  of  any  sub¬ 
structures  other  than  its  parent.  In  the  context  of  the 
flow-field  calculations,  the  air  mass  element  has  no  ac¬ 
cess  to  information  on  the  geometry  of  the  rotor.  This 
makes  it  virtually  impossible  to  make  any  assumptions 
that  would  allow  the  blade-element  contributions  to  the 
inflow  calculations  to  be  included  in  the  air  mass  ele¬ 
ment.  Any  assumptions  that  might  be  made  would  be  to 
the  detriment  of  the  generality  of  the  code.  Therefore, 
the  calculations  of  the  momentum  contributions  from 
the  actuator  disk  are  separated  from  the  blade-element 
calculations.  The  air  mass  element  represents  only  the 
flow-field  aerodynamics,  while  the  blade-element  aero¬ 
dynamics  are  isolated  in  the  aeroelastic  beam  element. 

THEORETICAL  DEVELOPMENT 

The  theoretical  development  of  the  inflow  equations  is 
dependent  on  three  components:  the  air  node,  the  air 
mass  element,  and  the  aeroelastic  beam  element.  The 
generalized  coordinates  that  are  used  by  GRASP  to  de¬ 
scribe  the  static  state  and  dynamic  perturbations  of  the 
induced  velocities  are  supplied  by  the  air  node.  The  air 
mass  element  performs  the  calculations  of  the  actuator- 
disk  contributions  to  the  inflow  equations,  while  the 
aeroelastic  beam  element  calculates  the  blade-element 
contributions. 

Under  the  assumptions  used  for  this  development,  there 
are  noncirculatory,  blade-element  contributions  to  the 
apparent-mass  terms  in  the  dynamic  inflow.  Some  re¬ 
cent,  but  as  yet  unpublished  work  indicates  that  the 
dynamic  inflow,  apparent-mass  terms  result  solely  from 
circulatory  effects.  If  this  can  be  verified,  some  of  the 
assumptions  used  in  this  analysis  would  have  to  be  re¬ 
vised. 


For  the  case  of  static  inflow,  generalized  coordinates  U f 
and  -yfr  are  used  to  represent  uniform  inflow  velocity 
and  radial  velocity  gradient  at  the  center  of  the  flow 
field.  The  other  two  coordinates  are  not  used.  Dynamic 
inflow  uses  only  generalized  coordinates  Uf,  qrfj,  and 
"lf3  to  represent  the  vertical  and  cyclic  velocity  pertur¬ 
bations. 

Air  Mass  Element 

The  air  mass  element  is  implemented  in  GRASP  to 
model  the  momentum  flow  of  air  through  the  disk  of 
a  helicopter  rotor.  In  this  element,  the  rotor  is  assumed 
to  be  an  actuator  disk,  and  the  flow  field  a  cylindrical 
region  surrounding  the  disk  (Fig.  1).  The  state  vector 
for  the  air  mass  element  is  made  up  of  the  generalized 
coordinates  for  a  single  air  node.  In  the  following  sub¬ 
sections  the  static  and  dynamic  inflow  models  developed 
for  the  air  mass  element  are  discussed. 

Static  Inflow.  In  the  static  case,  the  air  is  considered 
to  be  flowing  steadily  through  the  rotor  disk.  From 
momentum  theory  (Gessow  and  Myers,  1967),  the  dif¬ 
ferential  thrust  dT  acting  on  a  differential  annulus  of 
the  rotor  disk  is  related  to  the  induced  velocity  U  by 
the  equation 

dT  =  4xparU\V  +U\dr  (2) 

where  pa  is  the  air  density,  r  is  the  rotor  radial  coordi¬ 
nate,  and  V  is  axial  velocity  of  the  rotor  relative  to  still 
air  (positive  up).  The  total  virtual  work  done  by 
the  thrust  on  the  air  is 

«W=4*p0^  U6P\M  +  U\rdr  (3) 

where  e  is  the  root  cutout  radius,  R  is  the  rotor  ra¬ 
dius,  and  6P  is  the  virtual  displacement  of  the  air.  The 
expression  for  virtual  work  is  discretized  by  assuming 
that  the  induced  velocity  can  be  divided  into  a  uniform 
velocity  Of  and  a  radial  gradient  *tfT  so  that 


Air  Node 

The  induced  velocity  generalized  coordinates  are  intro¬ 
duced  into  GRASP  via  the  air  node.  These  generalized 
coordinates  are  defined  relative  to  an  inertial  frame  of 
reference  I,  and  they  define  the  inertial  air  velocity  at 
any  point  in  the  rotor  flow  field.  Given  that  is  an 
inertially  fixed  unit  vector  and  A  is  also  an  inertial  co¬ 
ordinate  system  with  its  origin  at  the  center  of  the  flow 
field  (Fig.  1),  the  induced  velocity  U^1  at  a  point  Q 


uQ‘  =  -(of  +  nt  +  S*  +  (i) 

where  r  is  the  flow-field  radial  coordinate,  and  R%f  is 
the  position  of  Q  relative  to  A  in  the  direction.  Uf, 
'll*,,  Tia,  &nd  Tij  are  the  air  node  generalized  coordi¬ 
nates. 


U  —  Of  +  *tfTr  (4) 

The  virtual  displacement  of  the  air  is  discretized  iden¬ 
tically.  Thus, 

6P  =  6Pf  +  6<frfrr  (5) 

When  these  expressions  are  substituted  into  the  expres¬ 
sion  for  the  virtual  work,  the  coefficient  of  6Pf  is  equal 
to  the  rotor  thrust  while  the  coefficient  of  64>\T  has  the 
units  of  moment,  but  no  real  physical  significance. 

Dynamic  Inflow.  The  model  for  the  inflow  dynamics  is 
taken  from  Pitt  and  Peters  (1981).  It  is  assumed  that 
the  freestream  velocity  of  the  rotor  relative  to  still  air 
is  spatially  and  temporally  uniform.  This  freestream 
velocity  is  augmented  within  the  cylindrical  region  of 


the  flow  field  by  the  steady-state  inflow  velocity  compo¬ 
nents  just  described.  Then,  infinitesimal  dynamic  per¬ 
turbations  to  the  inflow  are  induced  by  dynamic  per¬ 
turbations  of  the  rotor  thrust,  roll  moment,  and  pitch 
moment. 

The  virtual  work  for  the  unsteady  flow  of  air  through 
the  rotor  disk  is 


a:  2paU\V  +  U\6Pdil>dr 


US’- 


PaU6PdVe„ 


where  4>  is  the  rotor  azimuth  and  Ve(j  is  the  effective 
volume  of  the  cylindrical  flow  field.  This  statement  of 
virtual  work  produces  a  system  of  first-order  differential 
equations  may  be  converted  to  a  set  of  second-order 
equations  and  discretized  by  assuming 


U  =  Uy  +  +  P\  ~  ^i2rsin0  +  j>l3rcos\p  (7) 


=  -c,T“  (of  +  nf,  ■  K 

+*£<  +  «S<)  -  •*' 


The  relative  virtual  displacement  of  an  element  of  air 
with  respect  to  the  structure  SS<2„i  can  then  be  ob¬ 
tained  by  applying  Kirchhoff’s  kinetic  analogy  to  Eq.  9. 
All  C)  quantities  are  replied  with  6(  ),  and  all  velocity, 
angular  velocity,  and  velocity  gradient  symbols  are  re¬ 
placed  by  identically  labelled  virtual  displacement,  vir¬ 
tual  rotation,  and  virtual  displacment  gradient  symbols, 
respectively.  All  other  terms  are  then  discarded  from 
Eq.  9. 

The  magnitude  of  the  relative  wind  velocity  W  at  the 
aerodynamic  center  and  the  angle  of  attack  a  are  time- 
dependent  quantities  that  can  be  written  in  terms  of  the 
measure  numbers  of  the  relative  wind  velocity  vector. 
Since  this  theory  is  two-dimensional,  the  relative  wind 
velocity  and  angle  of  attack  depend  only  on  the  measure 
numbers  in  the  plane  of  the  blade  airfoil  cross  section; 
thus 


SP  =  &P*  -  6<t>i2rsintl>  +  60f3rcosV> 


=  /K',)j+K;)] 


•A 

where  Pl  is  the  vertical  perturbation  of  the  induced 

A  A 

velocity  at  the  center  of  flow,  <pl2  and  <pl3  are  the  cyclic 
perturbation  gradients  at  the  center  of  flow.  6P *  is  the 
vertical  virtual  displacement  of  the  air  at  the  center  of 
flow,  and  bd>\2  and  6<t>*3  are  the  cyclic  virtual  displace¬ 
ment  components  at  the  center  of  flow. 

Aeroelastic  Beam  Element 

The  aeroelastic  beam  element  is  the  primary  structural 
element  in  GRASP.  It  represents  a  slender  beam  that 
is  subject  to  elastic,  inertial,  gravitational,  and  aerody¬ 
namic  forces.  Hodges  (1985)  derives  the  elastic,  inertial, 
and  gravitational  forces  in  detail.  This  section  will  dis¬ 
cuss  the  derivation  of  the  aerodynamic  forces  as  they 
apply  to  the  induced  velocity  calculations. 

In  the  following  discussion,  the  symbol  Q  (for  quarter 
chord)  is  used  to  denote  the  aerodynamic  center.  The 
static  position  of  any  quantity  is  identified  as  (  )',  while 
(  )"  refers  to  the  instantaneous  position  of  the  dynamic 
motions  of  the  blade.  As  just  mentioned,  vectors  are 
denoted  by  the  underlined  symbol.  Measure  numbers 
of  vectors  associated  with  a  particular  set  of  unit  base 
vectors  are  subscripted  with  the  identifier(s)  for  that  set 
of  unit  base  vectors.  The  unit  base  vectors  used  in  the 
following  discussion  are  shown  in  Fig.  2. 

The  wind  velocity  vector  W®  at  the  aerodynamic  cen¬ 
ter  is  calculated  by  subtracting  the  inertial  structural 
velocity  at  Q"  (V^  f)  from  the  inertial  air  velocity  at 
Q"  ( U ^  ').  In  terms  of  the  inflow  generalized  coor¬ 
dinates  and  V 9  1 ,  the  relative  wind  velocity  measure 
numbers  associated  with  the  zero-lift-line  basis  vectors 


iy<3" 
_  "Z"! 


The  local  airflow  velocity  gradient  2  is  also  a  time- 
dependent  quantity  that  depends  on  the  relative  wind 
velocity.  The  subscripts  1  and  2  denote  the  gradient  in 
the  z3  direction  of  the  velocity  measure  number  in  the 
x i  direction.  This  velocity  gradient  can  be  shown  to  be 


qQ"  _  dwf, 

ki  Ttt  1  <1  —  />  rt 


which,  in  terms  of  the  inflow  generalized  coordinates 
and  V_Q  1 ,  is 


C£'.2  = 


-Z"A  dQ  * 
'22  nA2 


Ki  +  C33 


-*Z  A  nQ  ‘ 
'n  nA3 


dvM 


Like  the  virtual  displacement,  the  virtual  rotation  of  a 
structural  element  relative  to  the  air  6 T^„3  can  be  ob¬ 
tained  by  applying  Kirchhoff’s  kinetic  analogy  to  Eq.  13. 
This  is  accomplished  by  replacing  all  (‘)  quantities  with 
6  (  ),  replacing  all  velocity,  angular  velocity,  and  flow 


gradient  symbols  by  identically  labelled  virtual  displace¬ 
ment,  virtual  rotation,  and  virtual  displacement  gradi¬ 
ent  symbols,  respectively,  and  discarding  all  other  terms 
from  Eq.  13. 

Static  Inflow.  Since  the  relative  wind  velocity,  the  ve¬ 
locity  gradient,  and  the  angle  of  attack  are  all  time- 
dependent  quantities,  contributions  to  the  static  part 
of  the  virtual  work  are  obtained  by  separating  out  the 
static  terms  from  Eqs.  9,  13,  and  11.  The  static  magni¬ 
tude  of  the  relative  wind  is  then 


where  c  is  the  local  blade  chord,  W  is  the  magnitude 
of  the  relative  wind  velocity,  G®,, I2  is  the  flow  veloc¬ 
ity  gradient,  and  W®,,  t  is  the  flow  velocity  normal  to 
the  zero-lift  line  (Fig.  2).  The  lift,  drag,  and  moment 
coefficients  (e/,  ej,  and  cm),  respectively,  are  nonlinear 
functions  of  the  blade  angle  of  attack  a. 

Dynamic  Inflow.  After  the  static  quantities  in  Eqs.  9, 
11,  and  13,  have  been  removed,  only  the  dynamic  terms 
remain.  The  dynamic  part  of  the  magnitude  of  the  rel¬ 
ative  wind  is  written  as 


» = \j  (^z\y + a*) 


wQ'  i VQ"  a.  w®’  i VQ" 

llr  --  z  _ Z"2 

w 


(19) 


where  the  overbars  indicate  the  static  part.  Similarly,  ^  ^  dynamic  part  of  the  angie  of  attack  is 
the  static  value  of  the  angle  of  attack  is 


tan  a  = 


w Q' 

rr 


(15) 


I 1/Q"  _  iVQ" 
_  ’’  Z'l’’  Z"\  ”  Z'l"  Z"  1 

W2 


(20) 


The  expression  for  the  virtual  work  SM>  done  by  the 
aerodynamic  forces  over  the  length  of  the  beam  element 
is 


where  the  checks  indicate  the  dynamic  part.  Then,  the 
dynamic  terms  in  the  virtual  work  can  be  put  into  the 
form 


SW  =  [\-SS%:,Fz..,  +  Srf,3M)dx3  (16) 
Jo 

where  bS%f,  is  the  virtual  displacement  of  the  structure 
relative  to  the  air,  6T^..3  is  the  virtual  rotation  of  the 
structure  relative  to  the  air,  and  FZ"i  and  X  are  the 
applied  forces  and  moments  at  the  aerodynamic  center, 
respectively.  The  applied  aerodynamic  force  vector  F 
on  a  blade  section  (a  distributed  force  per  unit  length 
of  blade)  is  assumed  to  be 

F  =  +  Pb3"  +  Cncbf  (17) 


where  Cc  is  the  circulatory  lift,  V  is  the  drag,  and  £„e 
is  the  noncirculatory  lift. 

The  equations  that  define  the  aerodynamic  force  compo¬ 
nents  act  on  the  aeroelastic  beam  element  at  Q  and  are 
determined  from  a  quasi-steady  adaptation  of  Green¬ 
berg’s  thin-airfoil  theory  (Greenberg,  1947). 

Cc  =  ~PaW3CCi  +  ~PaCiWG<z„l2 
V  =  \p.W2ccd 
X  =  ^P«»Wcm 

-  ~p^  (wc%n2 + w?:x  + 

Cn'  =  \p*<?  (*«::,  + 

(18) 


where  and  <j5  represent  the  displavement  and  ve¬ 
locity  perturbations  of  all  of  the  structural  generalized 
coordinates.  From  this  expression  for  the  virtual  work, 
the  aerodynamic  contributions  to  the  aeroelastic  beam 
element  mass,  damping,  and  stiffness  matrices,  M,  C, 
and  K,  respectively,  can  be  determined  to  be 
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M  =  AFH 

C  =  AEH  +  AFG  +  BH  (22) 

K  =  AEG  +  BG+D 

Here  M  turns  out  to  be  symmetric,  but  neither  C  nor 
K  are.  Explicit  expressions  for  the  elements  of  M ,  C, 
and  K  can  obviously  be  obtained  by  substitution.  Such 
expressions  are  quite  long  and  complicated:  however,  in 
view  of  GRASP’s  method  of  evaluation  of  these  matrices 
numerically  from  Gauss-Legendre  quadrature,  it  is  not 
necessary  to  obtain  them. 

CONCLUDING  REMARKS 

The  method  used  in  GRASP  to  model  the  rotor  flow 
field  for  a  helicopter  has  been  described.  The  primary 
feature  of  this  implementation  that  differentiates  it  with 
other  approaches  is  the  separation  of  the  blade-element 
calculations  from  the  actuator-disk  calculations.  Also, 
this  method  incorporates  the  inflow  generalized  coordi¬ 
nates  in  the  state  vector  for  the  steady-state  problem, 
which  guarantees  full  coupling  with  the  structural  de¬ 
formations. 

Because  of  the  approach  used  to  implement  the  inflow 
calculations  in  GRASP,  it  is  also  possible  to  use  im¬ 
proved  flow-field  models  without  having  to  develop  an 
entirely  new  blade  element.  The  analyst  would  then 
have  at  his  disposal  a  prescribed-  or  free-wake  flow  field 
representation  as  in  Johnson  (1980),  or  perhaps  an  un¬ 
steady  flow-field  model  like  that  recently  developed  by 
Peters  and  He  (1987).  However,  since  the  current  ver¬ 
sion  of  GRASP  does  not  have  an  air  node  that  is  general 
enough  to  accomodate  the  different  sets  of  generalized 
coordinates  that  would  be  required  for  these  flow-field 
models,  an  improved,  generalized  air  node  would  need 
to  be  developed. 
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